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«U , Abstract 

"^ ' We perform a model dependent fit to recent data on charmless hadronic 



B decays and determine 7, the phase of V*f^. We find 7 = 114^21 degrees, 

C^ ' which disfavors the often quoted 7 ~ 60° at the two standard deviation level. 

We also fit for the form factors Fq^ and Aq^, and the strange-quark mass. 

They are consistent with theoretical expectations, although m-s is somewhat 

low. Such agreement and the good x^ for the fit may be interpreted as a 

confirmation of the adequacy of our model assumptions. 
PACS numbers: 



Typeset using REVT^ 



The measurement of a surprisingly large e' je value in 1999 is an exasperating reminder of 
how little we really know about CP violation in Nature. Within the Standard Model (SM) 
with 3 quark generations, however, there is a unique phase 0] in the Kobayashi-Maskawa 
(KM) matrix V , often defined as 7 = arg(\4*^) in the usual phase convention 0. At present, 
there is no evidence that this phase fails to account for CP violation phenomena. 

Two B factories, built to study CP violation in the B system, have just been completed. 
By comparing the time dependence of tagged B^ vs. B -^ J/ipKs decays, one can cleanly 
measure the CP phase in B^-B mixing, which, in the SM, gives sin 2/? where /3 = arg(V^^). 
Together with the demonstrated capabilities of collider detectors at the Tevatron, a precise 
measurement of sin 2/3 is assured within a year or two. The unitarity phase a can be 
measured via tt~^tt^ or tt+tt^tt^ modes but now appears to be more challenging because 
the TT+vr^ rate is smaller than expected, which in turn implies larger "penguin pollution". 
However, it is the phase 7 that is usually viewed as the most difficult. All suggestions so 
far require very high statistics or various technical challenges. In this Letter we exploit the 
emerging rare B decay data from CLEO and perform a fit 0] that, though model dependent, 
allows extraction of 7 with just 10^ B mesons using only CP-averaged rates. The goodness 
of fit and reasonableness of other fit parameters serve as checks on the adequacy of our 
model assumptions. Since neither vertexing nor tagging is required, this method will benefit 
from the improved statistics soon available from the CLEO upgrade as well. 

The measurement of sin 2(3 is often compared to a double-slit interference experiment, the 
two slits being B^ and B -^ J/ipKs decays. Charmless rare B decay rates, even when CP 
averaged, can also be viewed as double-slit experiments that in principle probe the phase 7. 
The present observed pattern that B — > Kit rates are larger than vrvr but comparable to 
p°7r^, p^tt''^ and ujtt^ implies that both tree (T) and penguin (P) amplitudes contribute to 
these rates, hence the double-slit analogy. Unfortunately, hadronization effects such as final 
state interactions (FSl) could dilute such interference effects. 

The Fleischer- Mannel bound Q on 7 is no longer effective since one now has R = 
T{B- -^ K-7i+)/T(B^ -^ K^Ti+) = LO ± 0.3, where T denotes the average of B and B 



widths. A more promising method is based on i?^, = r{B^ -^ K 'k^)/2V{B^ -^ K^n^), 
with some reference to tttt for control of model dependence. But with i?* = 0.75 ± 0.30 at 
present, one cannot set a useful bound on 7. More than an order of magnitude increase in 
data is needed for a restrictive measurement. In this Letter we take a more global view and 
perform fits which trade model independence for exhaustive use of available data. We also 
use only CP-averaged rates, since there is no sign of significant CP asymmetries ||^ and the 
errors are large. Asymmetries in addition are more sensitive to FSI than averaged rates. 

We shall assume that naive factorization {Nc = 3) is a good approximation, and use 
effective-theory matrix elements cross-checked by two groups 0,§|, ignoring annihilation 
type diagrams. We make a x^ fit of data to 7 and four other parameters. Factorization in 
two body rare B decays may be heuristically justified by the large energy release [^: final 
state mesons move away from each other so fast that they do not interact. Recent theoretical 
work suggests that factorization may be derivable from QCD in certain limits |]lOl- In our 



view, factorization provides a simple framework to describe hadronic B decays that is rich 
in predictions with a limited set of free parameters. It is therefore reasonable to use this 
framework when attempting a first global fit to the large number of results on charmless 



hadronic B decays now available. This work is motivated by Ref. [0, which pointed out 
that recent CLEO rare B data supports factorization if cos 7 < is taken. Indirect fits to 
the unitarity triangle find a 95% C.L. range for 7 of 44° - 75° III, 44° - 93° H, 41° - 



97° [0, and 36° - 97° fl^, depending in part on how conservatively the theoretical errors 



are handled. 

Let us illustrate the parameters that enter with B -^ K~7i~^, which is a 6 ^ suu 
transition under factorization. Ignoring annihilation terms, one has [0,H 

AK-n+ oc /kFq^'' (m| - ml) {V*^Vubai - V^lVtb [04 + aw + (ae + a8)Rsu]} ■ (1) 

We are free to fix |Vjs| = |V"cb| = 0.0381 since any uncertainty can be absorbed in form 
factors. The two relevant fundamental parameters are therefore iV^Mfe/Kbl and 7, and the 
latter clearly controls the interference between tree and penguin terms. The parameters 



Oj are related to short distance Wilson coefficients (WC) and evaluated within a QCD 
framework. They also depend on the scale parameter fif where factorization is operative. 
The values of a^ in the literature are still evolving as issues of scale, scheme and gauge 
dependence are addressed. We use two sets of a^ from Refs. ^ (AKL) and (CCTY). The 
dominant strong penguin coefficients are 04 and Og (—0.04 to —0.06), while the dominant 
electroweak penguin coefficient is og ~ —0.009 coming from the Z penguin. We use the Oj 
for b ^ s since the difference ioT b ^ d is small. 
In Eq. (|l|) one also has the factor 

Rsu = 2m^-/(mfe - niy){nis + m„) . (2) 

A similar factor Rsd-, which enters A^Ot^-, is taken to be equal to Rsu- This factor can be 
better understood as a product of two pieces: the factor l/(r7ife — ?7i„) = l/r/ifc balances 
against ni?g — m\\ and m\-/{ms + niu) is nothing but the nonperturbative part of the 
pseudo-Goldstone boson mass formula, which is well defined within QCD but not yet very 
well determined. Although Rg^ is technically related to an TTig-independent hadronic matrix 
element, in the form of Eq. (|^), it becomes a sensitive probe of iris in a way that is analogous 
to K ^ TiTi decay and e' /e. 

The factors Jk and F^'^ arise from evaluating hadronic matrix elements of four quark 
operators under factorization: The former comes from forming K^ out of the vacuum via 
the S7^75M current, the latter arises from the transition B^ -^ 7r+ via the ^7^6 current. 
While form factors are well defined, it is the reliance on models that causes us to lose track 
of the factorization scale /i/. Popular form factor models are the BSW model and light-cone 



sum rule (LCSR) evaluations. A recent compilation of models can be found in |jT6|, but we 
shall treat form factors as fit parameters. 

The criteria for choosing the decay modes to include in the fit are as follows. First, a 
central value branching ratio (BR) with statistical and systematic errors must be available. 
Second, we exclude VV modes such as uK* since there is insufficient data to constrain the 
extra form factors that enter. Third, we require that the experimental sensitivity (a few 



times 10~^ at present) to be comparable to the range of factorization predictions. Only the 
ujTT^ final state, with a predicted BR below 10~^, is removed with this criterion. Since this 
and other suppressed decays such as p^ir^, tt^tt^, 4>tt, KK and K*K may well be dominated 
by FSI from other charmless final states, factorization is less likely to work well. We therefore 
propose to exclude these modes even when suitable measurements become available. The 
only exceptions to these rules are final states involving r] and rj' . We prefer to apply the fit 
to predict their BRs rather than use them in the fit because the qq content and other issues 
of these mesons have recently been questioned. We note that the newly measured [|l^ rjK* 
modes, like rj'K, are larger than previous theoretical predictions 0,§]. 

We give the 14 measured BRs (averaged over B and B) that enter our fit in Table I, where 
we also give the fitted output. To limit the number of fit parameters, we use approximate 
relations as follows. We assume KM unitarity hence —V^gVth — |Kfe| and —V^^Vth = V*^Vcb + 
K*dKfe = -\Vcb\{X- e~''^\Vub/Vcb\), where A = |Ks|- Since A - \Vub/Vcb\ C0S7 > always, as 
noticed in Ref. [|1T|, T-P interference is opposite in sign for P-dominated and T-dominated 
modes such as K^n^ and 7r^7r+, leading to enhanced K^tc^'^ and suppressed 7r^7r+ for 
cos 7 < 0, in better agreement with data. The chiral relation m^./rn^- ^ {ms + mu)/{m(i + 
rriu) and the fact that m^ ^ m^^^ give Rsu — Rsd — Rdu, while Qij = —Rij for VP modes 
such as pvr, ujtt and uK. We use form factors at g^ = and Ff^ = F^^; Fq^/Fq"" = 1.13 
which is consistent with both BSW and LCSR models; A^^ = A^^; and Af^'* = 1.26 A^'' 
(used for predictions only). Surveying the amplitudes for modes in Table I, we find that just 
five parameters suffice for the fit: 7, iKife/Kbl, Rsu, F^^ and Aq^. 

The function minimized by the fit is 

X' = E((BRmeas " ^Kred) / <ej' + ((0-08 - |K./V;,|prcd)/0.02)2 , (3) 

i 

where we sum over the modes in Table I. The predicted BRs are calculated from formulas 
like Eq. (|I]) taken from Refs. 0,1^. We have checked that we confirm the Nc = 3 results 
of AKL and CCTY with the same input parameters. We take into account the full (asym- 
metric) experimental errors and correlations in K^n^ /ti^ti^, K^ti^/ti^ti^ and uK^ /ujtt~ 



measurements, where the correlation coefficients are —0.15, —0.29 and —0.17, respectively. 
The fit is able to nearly optimally use the information for each of these modes individually, 
though K/ti separation improvements in the next round of experiments will help in this re- 
gard. For simplicity, we assume that systematic errors have the same correlation coefficient 
as statistical errors, i.e. we apply the correlation coefficient to the total error with all errors 
combined in quadrature. If the best fit value is below the experimental central value, the 
low-side experimental error is used, and conversely the high-side. 

To understand the behavior of the x^ function in the 5D fit parameter space, its depen- 
dence on various fixed parameters or the exclusion of certain experimental measurements, we 
have explored many variants of our nominal fit. In all cases we find 7 > 100°. Our nominal 
fit results, with CCTY Oj values, are given in Table II. The x^ per degree of freedom (DOF) 
in the last column indicates the good quality of the fit. We choose CCTY rather than AKL 
as nominal only because of their claim of improved gauge dependence of the a^. The fit 
values for AKL (see Table II) differ only slightly from our nominal, and mostly because of 
the larger |a4e| found by these authors. We note that Rgu for AKL input should be smaller 
than the CCTY case since quark masses are defined at ^ = 2.5 GeV rather than rrih. We 
have also checked that the strong phases of 04^6 have little impact on our fitted 7 value. 

The x^ vs. 7 curves for the nominal fit with CCTY input are shown in Fig. [^. We 
note that our 7 value has a two- fold ambiguity since the fit is sensitive to cos 7 rather 
than 7. From the contributions from individual modes given in Fig. [l|(b), we see that the 
main discriminator for favoring large 7 comes from K~7t and tt^tt modes, and, somewhat 
surprisingly, the ujK~ and (j)K~ modes. The situation for (f)K~ is a result of the procedure 
of minimizing x^ for each 7 value. This induces an apparent sensitivity, due to changes in 
the other parameters, where there is no direct dependence. 

The error on \Vub/Vcb\ returned by the fit is only marginally better than the conservative 
range iKib/Kfel = 0.08±0.02 used as an additional term in Eq. (^. Sensitivity to | Kib/Kfel 
largely comes from B -^ p^Tf^, p^'K^ and cuvr^, all of which depend on Aq^ . Removing the 
constraint on | Kit/ Vet | from the fit gives higher iKife/Kbl with large errors and strongly 
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correlated with Aq'^ (and between Rgu and F^^'^ as well) but 7 = I21I24 degrees remains 
close to nominal (see Table II). 

The fit favors large Rsu (Rdu) since it facilitates the enhancement (suppression) of K^7r+'° 
(7r^7r+) modes. Furthermore, under factorization, the BRs for ujK modes are enhanced only 
for large Rsu such that 04 and ag penguin contributions do not cancel fully. We have 
checked that when the uK mode is removed from the fit, there is no significant change 
in 7, though Rsu drops to 1.69. Our nominal Rsu fit value implies rus = SSt^f, 67li3 
MeV at mb and 2 GeV scale, respectively. This is lower than what is commonly used 
in most previous calculations, but consistent with recent unquenched lattice results which 
give ms{2 GeV) = 84 ± 7 MeV |18|. In addition recent experimental results for e'/e can 
be reconciled better with theoretical predictions if a smaller value of rris is used [|l^. For 
comparison, the result for fixing Rsu = 1.21 {ms{mi)) = 90 MeV) is given in Table II. Note 
that Rsu is anti-correlated with og since only the product appears in the amplitudes. As for 
the form factors, our fitted F^'^ {Aq ^) is lower (higher) than but consistent with the LCSR 
result of 0.305 ± 0.046 (0.372 ± 0.074). 

Predictions from our fit for some selected modes are given in Table III. The agreement 
with the newly measured [0] rjK modes are rather striking. An enhancement factor of 1.7 
comes from Aq^' ^ 0.60 compared to LCSR value of Aq^' ~ 0.47, the rest coming from 
our low TUs. The rj'K modes are comparable in size to the observed rjK modes. Since 
we can account for less than half the rate of rj'K modes and the missing contribution may 
well be specific to the rj' decay modes, our predicted rj'K rates should be viewed with some 
caution. The p~7r° mode is suppressed by cos 7 < 0, smaller F^'" (which also suppresses 
K 7T modes) plus destructive interference between two terms because of low rris. The pK 
modes are enhanced by the low value of nis and larger Aq'' form factor, except for p^K~ ^ 
which is suppressed compared to p^K' by destructive interference between the strong ag 
and electroweak ag penguin terms (a similar effect suppresses K 7r° modes). 

As an aside, we give the "penguin pollution" as determined from our fit. Defining T (P) 



as the amplitude arising from ai_2 (fls-io); we find the ratio \P/T\ in ti^-k {p^n^) to be 
0.37 ± 0.04 (0.20 ± 0.04) for our nominal fit. For comparison, the CCTY result for A^^ = 3, 
7 ~ 65°, \Vub/Vcb\ = 0.090, ms{mb) = 90 MeV using LCSR form factors gives 0.20 (0.10). 

How do we reconcile with the usual fit to B and K data other than charmless rare B 
decays, which give a 95% CL. range that excludes cos 7 < [0? The removal of the second 



quadrant in these fits results mostly from combining recent bounds on Bg mixing from LEP, 



CDF and SLD, with lattice QCD results that relate Bd and Bg mixing parameters [|13 
Thus, our results suggest that Bs mixing could be very close to the present limit. 

It should be stressed that the goodness of our fit (see Table II) suggests that corrections to 
factorization may be small compared with the present experimental precision. It is reassuring 
that the hadronic parameters from our fit are not at variance with theoretical expectations. 
We note that 7 is the most stable parameter in the fit, with cos 7 < for all variations 
we have considered. This is because 7 directly controls the "double-slit" interference, while 
other parameters enter indirectly. We note that our larger value of 7 tends to reduce the 
value of sin 2/3. 

In conclusion, we have made a model-dependent determination of 7 = 114l|f degrees. It 
will be interesting to see if future, more precise measurements will confirm this result and 
the predictions in our tables. 

This work is supported in part by grants from the US DOE and NSC of Taiwan, R.O.C. 
We thank our CLEO colleagues for the excellent data. 



[1] 

[2] 

[s: 

[4] 
[5] 

[6] 
[7] 
[8] 
[9] 
[10 

[11 



[12: 

[13: 



REFERENCES 

M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 652 (1973). 

Particle Data Group, C Caso et ai, Eur. Phys. J. C 3, 1 (1998). 

Preliminary results are given in CLEO Collaboration, Y. Kwon et ai, CLEO-CONF 
99-14 and |hep-ex/990803U . 



R. Fleischer and T. Manuel, Phys. Rev. D. 57, 2752 (1998). 

M. Neubert and J. Rosner, Phys. Lett. B441, 403 (1998); M. Neubert, JHEP 9902, 14 
(1999). 



CLEO Collaboration, T.E. Coan et ai, CLEO-CONF 99-16 and |hep-ex/990802q , 



A. Ali, G. Kramer and CD. Lli, Phys. Rev. D. 58, 094009 (1998). 

Y.H. Chen, H.Y. Cheng, B. Tseng and K.C. Yang, Phys. Rev. D. 60, 094014 (1999). 

J.D. Bjorken, Nucl. Phys. (Proc. Suppl.) Bll, 325 (1989). 

T.W. Yeh and H.n. Li, Phys. Rev. D. 56, 1615 (1997); M. Beneke, G. Buchalla, M. Neu- 
bert and C.T. Sachrajda, Phys. Rev. Lett. 83, 1914 (1999). 

X.G. He, W.S. Hou and K.C. Yang, Phys. Rev. Lett. 83, 1100 (1999). The large 7 
possibility was first stressed by N.G. Deshpande, X.G. He, W.S. Hou and S. Pakvasa 
Phys. Rev. Lett. 82, 2240 (1999). 



F. Parodi, P. Roudeau and A. Stocchi, |hep-ph/9802289| and |hep-ex/ 9903063 



S. Mele, Phys. Rev. D 59, 113011 (1999). 



[14] S. Plaszczynski, in Proceedings of the Conference on HF8, Southampton, 1999 (to be 
published) . 

[15] A. Ali, D. London, Eur. Phys. J. C9 (1999) 687. 



[16] A. Deandrea, R. Gatto, G. NarduUi and A.D. Polosa, |hep-ph/9907225| ; and references 
therein. 



[17] CLEO Collaboration, S.J. Ricliichi et al, CLEO-CONF 99-12 and |hep-ex/ 990801^ 
M. Bishai et al, CLEO-CONF 99-13 and|hep-ex/99U80l8 



[18] S. Aoki, in Proceedings of the Symposium on LP99, Stanford, 1999 (to be published). 



10 







TABLES 








Table I. 


Measured BR (in lO'^ 


units) l^.l^ entering the fit and output BR from 


our fit, 


and for 


comparison (in parenthesis) the Nc = 3 values from Ref. [8| 


for 7 = 65° and LCSR form factors. 


Mode 


Measured BR 


BR from Fit 


Mode 


Measured BR 


BR from Fit 


K-TT+ 


18.8+1^ ±0.7 ±1.1 


20.9 (12.9) 


7r-7r+ 


4.7l||±0.5±0.3 


4.4 


(10.7) 


K-vrO 


12.lt|0t||±0.8 


11.8 ( 8.8) 


vr-vrO 


5.4+|J±1.5±0.3 


3.4 


( 5.7) 


K%- 


18.2^4:0 ±0.9 ±1.3 


19.5 (15.7) 


K vr" 


U.8tlitli ± 2.5 


8.0 


( 5.6) 


pOyr- 


15 ±5 ±4 


14.0 ( 8.1) 


pT^i 


35liJ ± 5 


39.0 


(41.3) 


UTr~ 


11.3l|^ ± 1.0 


11.5 ( 8.2) 


K*-Ti+ 


99+8+4 


5.4 


( 4.3) 


ujK^ 


3.il?:^±o.8 


3.7 ( 1.4) 


(pK- 


1.6+11 ±0.2 


3.0 


( 5.0) 


17^ 


10.0+^;^ ± 1.5 


4.3 ( 0.4) 


cfK' 


10.7l^;? ± 1.1 


2.8 


( 4.6) 



Table II. Results for fitted parameters. CCTY (AKL) implies use of a^'s from Ref. [§ (Ref. 0). 



Input Qi 


7 (degrees) 


\Vub/Vcb\ 


-tisu 


T^Bw 
^0 


A^'^ 


xiJBOF 


CCTY 


114^25 


0.087 ±0.016 


1 SQ+0-43 


0.26 ±0.04 


OAStlil 


11.3/10 


AKL 


105lg 


0.092 ±0.016 


1 7P+0.34 


0.23 ±0.03 


0.451°:^^ 


12.2/10 


CCTY 


124t| 


0.076 ±0.014 


1.21 (fixed) 


0.32 ±0.03 


0.58l°i? 


15.0/11 


CCTY 


1211^1 


0.105tHi (free) 


9 91+1.30 
^•^^-0.61 


p. 90+0.O6 

^•^"-•-0.07 


o.48t°:}g 


11.0/9 
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Table III. Predicted BR (in 10 ^ units) with CCTY OjS for some modes not used in the fit, and (in 
parenthesis) from CCTY for N^ = 3, j = 65°, iKfe/Kbl = 0.090, m,(m;,) = 90 MeV using LCSR 
form factors. For r/', r/ modes, we have used /"/, f5 = 64, 141 MeV, /^, f^ = 78, —113 MeV, and 
Fq"^', Fq"^ = 0.108, 0.121, respectively |J|. The /% effects are small. 



Mode 


Measured BR BRp^ed 


Mode 


BRpred 


Mode 


BRpred 


Mode 


BRpred 


rj'K- 


mtf ± 8 


32.4 (21.9) 


r^K- 


1.9 (1.7) 


p-.' 


6.5 (15.0) 


p+R- 


6.1 (2.0) 


r]'K 


88tJ^ ± 9 


28.1 (21.7) 


rj-r 


1.9 (1.0) 


K*-Ti^ 


5.1 ( 3.3) 


p^K- 


1.8 (0.5) 


T]K*- 


27.3l|^ ± 5.0 


16.2 ( 3.9) 


7]'K*- 


14.4 (2.0) 


K*^TT- 


3.5 ( 5.5) 


p K 


8.0 (0.4) 


rjK 


13.8^4:4 ± 1.7 


13.2 ( 4.3) 


rj K 


13.8 (1.0) 


^*0^0 


0.6 ( 1.3) 


p^K 


6.2 (1.1) 
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FIG. 1. (a) x^ — Xmin ^^- 7 from nominal fit, where Xmin = ^^-^ for 10 DOF; (b) Major x^ 
contributions from individual measurements. 
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